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PREFACE 


Hypersonic  strong-interaction  theory  provides  a  framework  for 
the  analysis  of  the  flow  past  a  slender  body  during  portions  of  its 
hypersonic  flight  through  the  atmosphere.  Performance  predictions 
and  scaling  laws  following  from  such  an  analysis  can  be  used  in  tra¬ 
jectory  and  discrimination  studies. 

In  this  Memorandum,  a  theoretical  study  is  made  of  the  laminar 
flow  past  a  very  slender  synmetrlc  body  when  the  Mach  and  Reynolds 
numbers  are  large,  in  which  case  Navler-Stokes  equations  are  reasonable 
models  for  the  equations  of  motion  for  a  viscous  heat-conducting  gas. 

It  is  well  known  that  the  flow-field  structure  for  two-dimensional 
bodies  exhibltr  Important  qualitative  differences  from  that  for  axi- 
symmetric  bodies,  and  one  objective  of  this  study  is  to  investigate 
these  differences  and  to  develop  a  consistent  theoretical  solution  of 
the  equations  of  motion.  Matched  asymptotic  expansion  techniques  are 
employed  to  develop  the  significant  features  of  the  flow  field,  and 
simple  expressions  are  presented  for  drag  and  for  heat  transfer  in  the 
extreme  case  where  the  body  diameter  is  much  less  than  the  thickness 
of  the  viscous  layer  supported  by  the  body. 

In  addition  to  its  relevance  for  discrimination  and  performance 
studies,  the  present  Memorandum  should  be  of  interest  to  those  Involved 
in  the  numerical  computation  of  viscous  flows,  since  it  provides  a  dif¬ 
ficult  test  case  that  can  be  used  to  make  needed  checks  on  numerical  work. 

This  study,  sponsored  by  the  Advanced  Research  Projects  Agency, 
is  part  of  continuing  RAND  Corporation  resesrch  on  hypersonic  gasdynsmics 
for  spplicstion  to  esrly  reentry  phenomenology. 


T 


SUMMARY 

The  Navier-Stokes  hypersonic  strong- inter act Ion  theory  is  presented 
for  the  flow  of  a  viscous,  heat-conducting,  compressible  fluid  past  a 
very  slender  axisymmetric  body,  where  D^,  the  measure  of  the  ratio  of 
the  radius  of  the  body  to  the  thickness  of  the  viscous  layer (s)  (which 
the  body  produces  and  supports),  goes  to  zero.  It  is  assumed  that  the 
fluid  is  a  perfect  gas  having  constant  specific  heats,  a  constant 
Frandtl  number,  o,  whose  numerical  value  is  of  order  unity,  and  vis¬ 
cosity  coefficients  varying  as  a  power,  u>,  of  the  absolute  temperature. 
Limiting  forms  of  the  solutions  are  studied  as  the  free-stream  Mach 
number,  M,  the  free-stream  Reynolds  number  based  on  the  body's  axial 
length,  R^,  and  the  modified  interaction  parameter,  » 

|M2 (*-“>) /R^og  (1/Db )}l/2,  go  to  infinity. 

By  means  of  matched  asymptotic  expansions,  it  is  shown  that,  for 
(1  -  u>)  >  0,  four  distinct  layers  make  up  the  region  between  the  shock 
wave  and  the  body.  The  self- similar  leading  approximations  for  the 
behavior  of  the  flow  in  these  four  regions  are  analyzed,  and  it  is 
found  that  an  arbitrary  body  supports  3/4 -power  viscous  and  shock 
layers. 
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SYMBOLS 


C£*Cq 


Sb,tb 

Sk,t:k 


friction  drag  coefficient  in  Eq.  (47a) 
pressure  drag  coefficient  in  Eq.  (47b) 

friction,  heat-transfer  coefficients  in  Eqs.  (46b)  and  (46c) 

pressure  coefficient  in  Eq.  (46a) 

dimensional  specific  heats 

6b/6,  thickness  ratio  parameter 

k, /k  ,  nondimensional  heat-conduction  coefficient 

characteristic  longitudinal  body  length 
1/2 

u_/ (YP^/P*,)  ,  free-stream  Mach  number 

p^/p^,  nondimensional  pressure 
P<nu®L^®>  free- stream  Reynolds  number 

r^/L,  nondimensional  radial  distance,  measured  from  axis 
of  symmetry  of  body 

nondimensional  independent  variables,  BVL  in  Eq.  (26a) 
nondimensional  independent  variables,  PVL  in  Eq.  (15) 

Tj/Too,  nondimensional  temperature 

Uj/u^,  Vj/u^,  nondimensional  axial  and  radial  velocity 


components 

(  2m  V/2 

\M  /Kj/  »  rarefaction  parameter 

(2u)  v  1  /  2 

M  /R^xj  »  rarefaction  parameter 

x^/L,  nondimensional  axial  distance,  measured  from  vertex 
of  body 

expansion  parameter,  PVL 
thickness  ratio  of  "effective  body" 
thickness  ratio  of  actual  body 


The  basic  nomenclature  is  given  here.  Additional  quantities  are 
defined  as  they  are  Introduced. 
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(yl)/(ytl),  Newtonian  parameter 

c{ (l+c)/(l-«)}  M  ,  stagnation  temperature  parameter 
[e{(l+c)/(l-c)}]l+U)  V2  Tb/6A,  parameter,  PVL 

nondimenslonal  second  viscosity  coefficient 
nondimenslonal  first  viscosity  coefficient 
Pj/p.,  nondimenslonal  density 


M^Cp  / Prandtl  number 

^  _  i 

[log (1/D. )]  ,  modified  thickness  ratio  parameter 


T  /0  ,  nondimenslonal  wall  temperature 

^  ^  1/2 

{m2^/^  log  (l/D^jj  ,  modified  interaction  parameter 


viscosity  exponent 


Subscripts 

b  body  viscous  layer  (BVL) 

h  inviscid  shock  layer  (ISL) 

ht  ISL-TL  matching 

k  principal  viscous  layer  (PVL) 

kb  PVL-BVL  matching 

kt  PVL-TL  matching 

t  viscous  transition  layer  (TL) 

w  body  surface 

1  dimensional  quantity 

•  free-stream  dimensional  quantity 
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I .  INTRODUCTION 

The  purpose  of  this  Memorandum  is  to  enlarge  upon  the  existing 
hypersonic  strong-interaction  i.heory  (HSIT)  for  viscous,  compressible 
flow  past  a  slender  axisymmetric  body,  whose  cross-sectional  radius  is 
much  less  than  the  thickness  of  the  viscous  layer (s),  which  it  produces 
and  supports. 

The  existing  HSIT  for  such  a  flow  regime  was  first  formulated  by 
Stewartson  (1964)  for  a  slender  co&  using  modified  von  Mises  coordi- 

v 

nates.  [Though  Stewartson's  formulation  is  for  a  cone,  his  approach 
is  a  modification  of  the  formulation  for  a  "3/4-power"  body  --  the  body 
for  which  similarity  solutions  have  been  shown  to  exist  when  the  body 
radius  and  the  thickness  of  the  viscous  layer  are  of  the  same  magni¬ 
tude  (Yasuhara  1956,  1962;  Luniev  I960) . )  Succeeding  analyses  by 
Solomon  (1967)  and  Ellinwood  and  Mirels  (1967)  have  followed  the  path 
laid  down  by  Stewartson,  employing  the  Lees-Dorodnitsyn  "boundary  layer 
coordinates,"  which  are  not  especially  appropriate  for  the  flow  problem 
under  consideration  (as  the  latter  authors  point  out). 

In  the  Stewartson-Solomon-Ellinwood  and  Mirels  (SL'EM)  formulation 
a  shock  wave  is  postulated,  and  between  the  shock  wave  and  the  body 
there  are  three  principal  zones  of  interest:  (1)  the  thin  inviscid 
zone,  or  inviscid  shock  layer  (ISL),  described  by  hypersonic  small- 
disturbance  theory  (HSDT),  bounded  on  one  aide  by  the  shock  wave  and 
on  the  other  by  the  viscous  layers  on  the  actual  body;  (2)  the  thin 
outer  viscous  zone,  or  principal  viscous  layer  (PVL),  comparable  in 
thickness  to  the  ISL,  in  which  the  fluid  velocity  and  temperature  are 
slightly  changed  from  their  ISL  values;  and  (3)  the  inner  viscous  zone, 


or  body  viscous  layer  (BVL)  --  which  is  implicitly  of  the  thickness  of 
he  body,  end,  hence,  thinner  than  either  the  ISL  or  the  PVL  —  in  which 
the  sheer  end  heat -conduction  terms  dominate  the  tangential  momentum  and 
energy  equations,  respectively. 

In  the  SSEM  formulation,  the  matching  of  the  solutions  for  the 
various  zones,  or  layers,  is  not  treated  formally  in  either  the  ISL- 
PVL  or  the  PBL-BVL  case. 

In  the  following  sections,  a  rigorous  and  more  complete  treatment 
of  the  very  slender  axisymmetrlc  body  HSIT  is  presented  for  the  case 
where  self-similarity  exists  throughout  the  flow  field.  In  Section  II, 
the  Navier-Stokes  equations  are  given  in  cylindrical  polar  coordinates. 
In  Section  III,  the  similarity  HSDT  for  the  ISL  is  given  under  the  pre¬ 
sumption  that  the  "effective  body,"  made  up  of  the  viscous  layers  and 
the  actual  body,  is  a  power-law  body.  The  treatment  is  essentially 
that  given  by  Stewartson  (1964).  In  Section  IV,  the  similarity  formu¬ 
lation  of  the  "Oseen-like"  PVL  is  presented.  This  formulation  follows 
quite  closely  that  for  Stewartson's  outer  viscous  zone,  except  that  it 
Introduces  a  modified  radial  variable,  which  simplifies  the  analysis 
of  the  axisymmetrlc  Navier-Stokes  equations.  In  Section  V,  the  simi¬ 
larity  formulation  cf  the  "Stokes- like"  BVL  is  given.  As  in  the  case 
of  the  PVL,  the  proper  choice  of  the  modified  radial  variable  greatly 
simplifies  the  analysis  of  this  layer,  which  follows,  in  part,  that  of 
Stewartson's  inner  viscous  zone.  The  details  of  the  PVL-BVL  matching 
are  also  presented. 

It  is  determined  that  the  theory  presented  is  valid  for 
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Further,  it  is  found  that,  to  the  order  of  approximation  considered 
for  the  flow  quantities,  the  PVL-BVL  matching  requires  the  "effective 
body"  shape,  i.e.,  the  outer  edge  of  the  PVL,  to  be 

r  -  ^  “  Nh(c)*  const. 

but  puts  no  requirement  on  the  actual  body  shape,  which  may  be 
represented  as 


r  ■  (x),  Ffa(x)  "  arbitrary  fnc(x) 

The  pressure,  friction,  heat  transfer,  and  drag  coefficients  are 


where  C  ,  C,  .  C  -  fncs(«,  a,  tw,  T  /0_)  »  const.  These  expres- 
p,o*  f to '  q,o  w  5 

sions  for  the  coefficients  are  equivalent  to  those  found  through  the 
SSEM  formulation. 

In  tha  Appendix,  the  formulation  of  the  required  viscous  transition 
layer  (TL),  intermediate  to  the  1SL  and  PVL,  is  given  (cf.  Bush  1966). 
Tha  1SL-TL-PVL  matching  for  the  temperature  is  demonstrated. 

A  schematic  diagram  of  the  hyperaonic  strong-interaction  layers 
for  flow  paat  a  very  slender  axlaymmetrlc  body  is  given  in  Fig.  1. 


II.  THE  EQUATIONS  OF  MOTION 


Consider  the  flow  of  a  viscous,  compressible  gas  past  a  very 
slender  axisymmetric  body.  Let  ■  Lx  and  r^  ■  Lr  represent  the 
cylindrical  polar  coordinates  along  the  axis  of  symmetry  from  the 
vertex  of  the  body  and  normal  to  this  axis,  respectively.  The  length 
L  is  chosen  so  that  x  is  of  order  unity  in  the  region  where  the  strong* 
interaction  theory  is  valid. 

Under  this  formulation,  the  equation  of  the  surface  of  this 
slender  body  is 

r  -  $bFb(x)’  with  6b  “  6Db  <<:  6  <<:  l>  Fb(x)  "  0(1) 

where  6  represents  the  scaling  of  the  effective  thickness  of  the  vis¬ 
cous  layer (s)  supported  by  the  body.  The  velocity  components  in  the 
x,-  and  r, -directions  are  u,  ■  u  u  and  v,  ■  u  v,  and  the  pressure, 
temperature,  and  density,  respectively,  are  p^  ■  p^p,  T^  ■  T^T,  and 
Pj  ■  p^P,  where  ub,  p^,  T^,  and  p^  are,  respectively,  the  velocity  in 
the  x^-direction,  pressure,  temperature,  and  density  in  the  undis¬ 
turbed  region  upstream  of  the  body. 

A  perfect  gas  (p  -  pT)  is  assuaed,  having  (1)  constant  specific 

heats,  c  end  c  ,  with  v  ■  (c  /c  )  ■  const.;  (2)  a  constant  Prandtl 
V1  ^1  ^1  V1 

number  of  order  unity  (o  ■  const.  -  0(1));  and  (3)  its  first  and  second 
viscosity  coefficients  proportional  to  a  power,  us,  of  the  absolute 
temperature  (p^  -  -  pjr**’,  with  1/2  £  us  <  1;  X^  «  p^X  -  Kp^p  -  Kp^T10, 


K  -  const.  ■  0(1) ] . 
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The  Navier-Stokes  equations  of  motion  in  cylindrical  polar  co 
ordinates  for  the  flow  of  such  a  gas  are 


(pur)  +  I;  (Pvr )  -  0 


(1) 


P 


(• 


du  du\  ,  1  -  e  1  dp 
dx  v  dr/  1  +  e  Jl  dx 


(lU>{(2  +  K) 


(2) 


P 


(u 


dx 


dv\ 

dr/ 


+ 


1  -  e  1  dp 
1  +  e  M2  dr 


+ 


2 

1  +  e 


(4) 


where  e  -  (y-1)/(y+1)  ■  0(1),  the  Newtonian  approximation  (of  e  «  1) 
not  being  invoked;  -  (P^u^/yp.)  »  1;  ©<,  “  e{ (l+e)/(l-e)}M^  »  1; 
and  -  (p.u.L/yJ  »  1. 


n 


■)  H**X  *{T 
Qt 
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III.  THE  INVISCID  SHOCK  LAYER 


According  to  present  hypersonic  strong-interaction  theory  (HSIT), 
the  slender  body,  whose  surface  is  given  by  r  ■  6jjFjj(x),  combines  with 
the  thin,  viscous,  heat-conducting  layer(s)  at  the  body  surface  to 
disturb  the  uniform  external  flow.  This  combination  of  the  body  and 
the  viscous  layer(s),  whose  outer  edge  is  given  by  r  ■  6F^(x),  with  6, 
the  thickness  parameter  of  this  combination  (for  the  flow  regime  under 
consideration,  6  ■  6D,  «  6  «  1)  acts  as  a  slender  "effective 

D  D 

body,"  producing  an  oblique  shock  wave,  r  ■  6F^(x)  >  6F^(x),  and 
an  inviscid  shock  layer  (ISL)  between  the  shock  wave  and  the  "effec¬ 
tive  body."  For  this  ISL,  whose  lateral  extent  is 


6Fk(x)  *  r  s  6Fsh(x) 


the  approximations  of  hypersonic  small-disturbance  theory  (HSDT)  are 
assumed  to  apply. 

Introducing  the  coordinates  x^  and  r^,  defined  by 


r 

6 


(5) 


the  HSDT  expansions  of  the  dependent  variables  have  the  form 


l  +  l  uh  +  ... 


6v^  +  ... 


2  2 

MTph  +  ... 


(6) 


with  fh  -  fh(xh,rh) 
HSDT,  taking 


t  -  es6‘Th+... 

0(1).  Thus,  the  first-approximation  equations  of 


f 


A. 


1 


are 


The  Rankine-Hugoniot  shock  relations,  which  determine  the  boundary 
conditions  for  the  flow  quantities  of  Eq.  (6)  at  rh  -  pBh(*h)»  are 


(V.h  ■ 

-<1  * 

‘V-h  ■ 

<*  '  •>'*<%>  ) 

(ph>8h  • 

(1  + 

c)[F8'h(Xh)]2, 

<V.h  * 

<*  •  ‘)CFs'h<*h>^j 

To  obtain 

similar  solutions  to 

the  general 

HSIT  problem  for  a 

very  slender  axisymmetric  body,  it  is  necessary  to  obtain  similar 

2  2 

solutions  for  the  ISL.  For  H  6  »  1 ,  the  flow  in  the  ISL  is  self- 

similar  if  the  shock  associated  with  this  flow  is  described  by  a 
power  law;  i.e.,  if 


In  the  notation  of  Section  IV,  [  (c{  (l+«  )/(l-c)})^+Uk^U>/RT  ]  ■ 

2  1+U)  " 

A|t(6  /o^)  •  In  the  formulation  of  Section  IV,  it  is  necessary  that 

■  0(1),  (6^/a^)  «  1;  hence,  the  requirement  expressed  above  for  the 
ISL  is  consistent  with  the  requirements  of  the  PVL. 
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F.h(xh>  *  <9> 

For  such  similar  solutions,  the  appropriate  independent  and  dependent 
variables  are 

?h  ■  v  \  ■  <"» 

-h 

uh  ■  •<>  -  ‘)"25h2<1'">uh<\)’  \  ■  <■  -  ‘K<l-m,VV  \ 

;  do 

"h  ■  u  +  tk  ■  a  -  o-!!;2(1'\(\)) 


When  Eqs .  (10)  and  (11)  are  introduced,  the  equations  of  motion, 
Eq.  (7),  reduce  to 


when  Wh  -  [(1  -  t)V^  -  T^].  The  boundary  conditions  for  these 
equations,  obtained  from  the  shock  relations,  are 
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The  solutions  of  Eq .  (12),  subject  to  Eq.  (13),  determine 

Nh  “  Nh(«,  m)  <  thc  valuc  of  \  at  the  outer  edge  of  the  "effective 

body,"  where  ■  0  (cf.,  e.g.,  Lees  and  Kubota  1957).  For  1/2  <  m  <  1, 

as  Tl,  -  N.  ,  the  similarity  solutions  yield 
'h  h 


P.  -  P.  (e,m)  +  ... 
h  h,  o 


vh "  ' C)1  + 


(14«) 


-E 


h  -  elii0(..m)(vv 


E  ■  •  *] 

Hence,  the  expressions  for  ph,  vh,  and  T^,  as  rh  -  are 

.  .  2.,  -2(1 -m) 

ph  -  ( 1  +  C  )m  h  +  .  . . 

_  .  v  -2(1 -m) 

-  nh(c,m)5h  +  ... 


-(!-■)  + 

.  T  •  *  i 


-"Vt» 


-  (1  -  t)m\  t 


(lAb) 


■2(1 -m) 


-E 


‘hi'*  -'“  ''h.o^h  (\  ’  Nh)  + 

-  *h(«.-)«;2<1'B)(\  -  vB+  ••• 


For  m  -  3/A,  c  -  1/6,  ^(1/6,  3/A)  -  0.875,  and  iyi/6,  3/A)  -  0.5A9; 


E(l/6,  3/A)  -  0.313. 
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IV.  THE  PRINCIPAL  VISCOUS  LAYER 


Just  as  there  is  a  principal  viscous  layer,  usually  designated  as 
the  viscous  boundary  layer,  which  acts  as  an  "effective  body"  in  the 
HSIT  for  flow  past  a  flat  plate,  it  is  postulated  that  there  is  a  * 
principal  viscous  layer  (PVL),  having  a  thickness  ratio  of  0(6)  (with 
6.  ■  6D.  «  6  «  1),  which  acts  as  an  "effective  body"  in  the  HSIT  for 

D  D 

flow  past  a  very  slender  axisymmetric  body  with  a  thickness  ratio  of 

0(6.  ).  As  in  the  case  of  the  HSIT  viscous  boundary  layer  for  a  flat 
b 

plate,  the  HSIT  PVL  for  a  '»ery  slender  axisymmetric  body,  which  sup¬ 
ports  the  ISL,  is  considered  to  be  a  high  temperature,  low  density 
region,  across  which  the  pressure  is  constant,  at  whose  outer  edge, 
r  *  6Fr(x)  -  6N^xm,  the  flow  quantities  u,  T,  v,  and  p  have  the  be¬ 
havior:  u  «  1,  T/0g  -  0,  v  -  6mNhx  ^  m\  and  p  -  M^I^x 
where  and 

The  appropriate  coordinates  of  the  PVL  are  taken  to  be  s^  and 
tfc,  given  by 


are  constants  defined  in  Section  III. 


s 


k 


x, 


(15) 


From  the  definition  of  t.  ,  it  follows  that 

k’ 


tk  -  0  as  r  ->  6Fk; 


®  as  r  -  6D.F. 

b  b 


The  expansions  for  the  flow  quantities  in  this  layer,  following,  in 
part .  Stewartson  (1964),  are 
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u  -  1+Vk+ 


v  -  6F'  exp  (-tk)v.  +  ...  -  6(mNhsk(1'm))  exp  (-tk)vR  + 


T  ■  ®sak^k  +  *  * ‘ 


2  2 

MTpk  +  • 


with  <*k  *  parameter  much  less  than  0(1)  to  be  determined,  and  fk  * 

£k(W  ■  »<«)• 

For  these  representations,  with 


(•M)1"  5  «T  ■  «a>.  »2  «  «  •’ 

V 


the  leading  terms  in  the  equations  of  motion  for  the  PVL  are 


-  o,  i..„  P 


k  ■  W  -  Vk2<1'm) 


/'  •  V 

1  •  vk\ 

_  l 

l  Tk  ) 

2 

\  Tk  ) 

m 

(17a) 


a 

)  -  -  i  [2«"  • » t +  s-  ^  (t)|  <17b> 


*The  requirements  of  (1)  ■  0(1),  (2)  6  2  «  ar^  «  1, 

and  (3)  M262  »  1  yield  the  following  individual  inequalities: 

m-2/(1-ui)  <K  b  ^  l  m-A/(1-uj)  <<;  <<;  x 

k 

1/2 

1  «  M2u,«  «  M2(2+U)),  or  M-2  «  VL  -  (M2u,/Rl )  «  1 


Further,  the  modified  interaction  parameter,  x* >  obtained  from  (1)  and 
(3)  is 


X'  “  {m2(2+uj)oJ+("/rl} 


»  1 
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_  m2  4m-3  , 

"kVk  exp  <■ 


_  m2  4m-3  . 

nh»h»k  < 


sk  Suk  I1  '  vk\  Suk 

m  a  d  /_ u)  ^Uk\ 
k  atk  \  k  atk/ 

s.  ST.  /I  -  v.  \  aT.  . 

'2tk)k^  +  m(~)^;+i  +  « (1 

.  _a_  (jw 

a  atk  \  k  atj 


(17c) 


(17d) 


Equations  (17b-d)  must  satisfy  the  boundary  conditions 


ui,(sir»  0)  "  ^ir(8w»  m  0,  vl-(si,»  0)  ■  1 


at  the  outer  edge  of  the  PVL  in  order  that  the  PVL  solutions  may  match 
to  those  of  the  1SL. 


Uk(8k’tk)  "  Ak<8k)Uk(tk) 


Tk*8kftk*  "  M^^k^k* 


vk(8k'tk)  "  W 


(19a) 


where 


Ak<8k> 


2  l/(l+u>) 

/mnhNh  \  (4m-3)/ (1+tw) 

\  Ak  /  k 

,  2  x/ (l+w) 

m 


(19b) 
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l 


the  self-similar  ordinary  differential  equations  and  boundary  oonditions 
governing  the  flow  in  the  PVL  are 


(20c) 


Uk(0)  *  ®k(0)  "  °»  \(°)  “  1  (21) 


Indeed,  it  is  found  that,  in  conjunction  with  Eqs .  (20s)  and  (21), 
Eqs.  (20b)  and  (20c)  yield,  upon  Integration, 


dU, 

— -  -  2 
k  dt,  * 

k 


(V)/' 


\xp  (-  2~k)  ^  dtk 
k 


U, 


-  e 


(Hr)  l1  ■  <-  2tk>J +  ^  •  V  “■>  <•  V 

(22a) 


i0^ 
o  wk  dt, 
k 


^®k  T /2m  -  1  •  2e  /I  •  m  T, 

dt”  "  *  rn  {—)}  i1  ■  e*p  <-  2VJ 


+  (1  -  Vk)  exp  (-  2tk) 


(22b) 


Near  the  outer  edge  of  the  PVL,  where  tk  -»  0,  it  is  found  that,  for 
( 1  -  uu)  *  0(1)  >  0,  the  asymptotic  solutions  for  the  flow  quantities  in 
Eqs .  (20)  are 


r 
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e.  -  ®„  t2'<i*»>  +  „ 

k  k,o  k 


\  +  - 


V.  -  1  -  V.  t,  +  ... 

k  k,  o  k 


(23a) 


where 


ek,»  •  -Is  (H^X2<2"  -  »  -  *  rh-  <l  -  "X1  -  -»]} 


vsl/(l+U)) 


\t0  ■  (1  -  "0(1  -  {2(2m  -  1)  -  ~  3>}  -  (1  -  ®)(1  -  »)]  1  efcj 


v  .  I  (i±»)  r2(2m .  i) . 

k,o  m\l-u)/L  1  +  a)  J 


(23b) 


Wh^n  m  -  3/4, 


©, 


k,o 


A  comparison  of  Eqs.  (14)  and  (23)  shows  clearly  that  the  func¬ 
tional  behaviors  of  the  temperature  in  the  ISL  and  the  PVL,  as  -*  N^  # 

and  tfc  -♦  0,  respectively,  do  not  permit  direct  matching  between  these 

* 

two  layers,  and  the  introduction  of  an  intermediate  transition  layer 
is  necessary  to  complete  the  matching. 

Near  the  inner  edge  of  the  PVL,  where  t^  -♦  °°,  it  is  seen  from 
Eqs.  (22a-b),  that 


The  details  of  the  formulation  of  this  transition  layer,  which  is 
consistent  with  that  introduced  in  the  HSIT  for  flow  past  a  flat  plate 
(Bush  1966),  are  presented  in  the  Appendix. 
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0 


dU. 

U!  _ k 

k  dt, 


exp  (-  2tR)  dt 


e(l 


U) 


d© 


k  dt 


k 

k 


i  [(2"  -  »  +  TT7  t1  *  ">] 


Hence,  the  shear  and  heat -conduction  terms  are  found  to  dominate  the 
momentum  and  energy  equations,  respectively,  yielding,  in  this  limit, 
the  following  asymptotic  solutions  of  the  similarity  equations: 


0. 


■  *  -=>T  *  es°'k\ek.-ti/<1+u’)  +  - 


If  _l/(14y))  .  v  1  .  .  II  l/(l+u>)  , 

Vck  +  —=•>“  ■  1  +  Wk.A  +  — 


\  “  vk,-+— ‘  5('nNh,k<I"n))vk,~  exp  V  +  ••• (24a) 


where 


0, 


k|°° 


-  (i  +  '»))[( 2m  -  1)  +  y~7  (1  -  m)]} 


1/(1+U)) 


J  _  1  +  (1) 

k,®  m 


00  U,  -I 

p  f*  ir  _  (JJ 

[2(2m  -  1)  exp  (-  2tR)  g-  dtR  -  s(l  -  m)J  0R>o 


v  =  (-Li- 2 1)  +  — L_ 

k,®  \  m  /  1  +  ui  \  2m  / 


(24b) 


With  the  possible  exception  of  the  asymptotic  solution  for  v, 
the  solutions  of  Eqs .  (24a-b)  do  not  satisfy  the  usual  nonslip, 
temperature-specified  boundary  conditions  at  the  body  surface, 
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VI 

Db 

(25) 

Therefore,  the  analysis  of  an  additional  layer,  interior  to  the  PVL, 
in  which  t^  ■  o([log<l/Db)])  -  ®,  is  introduced  in  the  next  section 


u*v  -  °*  T  - 

s 


specified  fnc(s^), 


as 


tk  -  log 


in  order  to  describe  the  adjustment  of  the  flow  in  the  vicinity  of 
the  body  surface  to  satisfy  the  surface  boundary  conditions  of  Eq .  (25). 
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V.  THE  BODY  VISCOUS  LAYER 


For  the  viscous  layer  adjacent  to  the  body  surface,  referred  to 
in  this  Memorandum  as  the  body  viscous  layer  (BVL) ,  the  appropriate 
coordinates  are  s^  and  t^,  given  by 


Tb 108 


'*Fb(x)l  ^ ,  /F b\‘ 

~  J  ■  Tb  [s + 108  (f^I 


(26a) 


where 


Tb  "  [log  (d:)]  ^  1 

D 


(26b) 


From  the  definition  of  t^,  it  follows  that 

Fb 

tb  -♦  1  as  r  -*  6DbFb;  tfa  -♦  Tb  log  ~  -♦  0  as  r  -♦  6Fj 

The  expansions  for  the  flow  quantities  in  this  layer  [cf. 
Stewartson,  et  al.  (1964)]  are 


u  ■  u,  +  . . . 


T  -  esTb  +  ... 

2  2 

p  ■  M  6  pb  +  ... 


(27) 


with  fb  -  f b (sfa ,  tfa)  =  0(1). 
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Subject  to  later  verification  as  to  consistency,  for  t 
6"Va^+UI,  and  (Qf^+UJ/T^)  exP  (“2/t^)  «  1,  the  leading  terms 
equations  of  motion  for  the  BVL  are 


3tb 


0, 


i.e. 


W 


nb»b 


The  boundary  conditions  to  be  satisfied  at  the  body  surface 


“b<vl>  ■  °-  Wl>  ■  \ 


Taking  the  wall  temperature  to  be  constant  so  that 


cp  =  @  a  const, 
w  b,w 


and  taking  the  flow  variables  to  be  of  the  form 


W  tb)  ■  W 
VV  V  "  eb(tb) 


vb(V  eb> 


w 


b’ 

in  the 

(28a) 


ub  (28b) 

(28c) 


(28d) 

are 

(29) 


(30a) 


(30b) 


I 
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the  self-similar  equations  and  boundary  conditions  governing  the  flow 
in  the  BVL  are 


V 


b 


U 


b 


ub(U  ■  0,  eb(i)  -  0b>w 

The  first  integral  of  Eq.  (31b)  is 


0 


dUK 

ui  _ b 

b  dt. 


-  S. 


(31a) 

(31b) 

(31c) 

(32) 

(33a) 


where  S^,  the  shear  function,  is  a  positive  constant.  The  first 
integral  of  Eq .  (3lc)  is 


0 


d0, 
a>  _ b 

b  dt. 


(33b) 


where  ^  ,  the  maximum  value  of  0^  in  the  BVL,  is  a  positive  constant. 

The  positive  branch  of  Eq.  (33b)  represents  the  temperature  field 

between  the  point  of  maximum  temperature,  t^  ■  t^  and  the  outer 

edge  of  the  BVL,  t^  -♦  0;  the  negative  branch  of  Eq.  (33b)  represents 

the  temperature  field  between  the  point  of  maximum  temperature, 

t.  *  t,  ,  and  the  body  surface,  t,  -»  1. 
b  b,m  b 
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For  -  1  and  0^  -  0  as  -•  0,  Eqs.  (33)  combine  to  yield  the 
Crocco-like  second  integral, 


U 


b 


(34) 


whose  negative  branch  is  valid  for  0  £  t.  <  t.  ,  and  whose  positive 

b  b  ,m 

branch  is  valid  for  t.  <  t.s  1  .  From  an  evaluation  of  the  positive 

b  ,m  b 

branch  at  t,  M  1,  i.e.,  for  U.  ■  0,  ®.  “  0.  ,  it  is  determined 

b  b  b  b,w 

that 


i  a 

4  1  +  e 


1  +  e 
a 


>  0 


(35) 


When  Eq .  (33b)  is  solved  to  determine  the  solutions 

UL(tL)  and  V, (t, )  follow  directly  from  Eqs.  (34)  and  (31a). 
b  b  b  b 

From  direct  integration  of  Eq .  (33b),  taking  ©b(0)  ■  0  and 
qb(1)  ■  0^  w,  it  is  found  that  the  shear  function,  S^,  is 


*To  give  an  idea  of  the  magnitude  of  Sb>  it  is  noted  that 


S 


for  a 


1 

6 


dZ 
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Thus  ,  near  the  outer  edge  of  the  BVL,  where  t^  -*  0,  it  is  seen 
from  Eq .  (33b^  that 


uu 

% 


fb„_°_s  C,  +  U_e  e  , 

dt,  1  +  e  b  \  a  b,vj 


and,  hence,  it  follows  that,  at  the  outer  edge  of  the  BVL, 

l/(l+a))  ,  v  T  -  ft  ©  i/d+O)) 

%  =  Qb,olb  +  ••’  *=>  T  6S°b,otb  ••• 

1/(1+id)  .  v  ,  i  d/(l+UJ)  . 

"b  ■  1  -  Ub,oV  +  •••  ->  u  1  •  Ub,o‘b 


V.  -  V.  „  +  .  • 

b  b  ,o  ' 


v  ■  6(1  -  m)  —  V  exp 

sb  b,0 


("  ’b)+  - 

■  5["Nh%'(1"n>]  [(H^K. 


o  I  exp  <-*k>  +  •• 


where 


P) 


b  ,o 


b.o 


> +  “)Sb(rTrXl  +  ^  eb,w)j 


-1/(1+U)) 


(1  +  uu)S  ci  a 
b  b  ,o 


(37a) 


V «  =  1 

b  ,o 


(37b) 


To  verify  the  similarity  formulations  of  the  PVL  and  BVL,  at  least 
insofar  as  these  two  layers  are  interrelated,  it  must  be  shown  that  the 
solutions  of  Eq .  (24a-b)  match  to  those  of  Eq .  (37a-b).  The  solutions 
of  the  PVL  and  the  BVL  are  matched  through  the  introduction  of  the 
intermediate  limit,  lim^,  defined  by 


r 


8.  ,  ■  x  f  ixed , 
kb 


'kb  ■  Tkb  108 


fixed ; 


Tb  <<:  Tkb  <<:  1*  Fkb^  “  arbitrary  fnc(x)  «  0(1) 


In  this  limit, 


F 

'k  ■  tTT  ['kb +  Tkb 108  Grr)]  -  tz  'kb  - 


'b  ■  r:  [\b +  \b 108  (ft)]  -  rr  'kb  - 0 


kb  “•  kb 


The  PVL-BVL  temperature  matching  requires  that 


,[0sVk(8kb-  +  •••)  +  •  •J  -  9STb (’kb -  +  •■•)  '  ••• 


(39a) 


From  Eqs .  (24)  and  (37),  it  is  seen  that  Eq .  (39a)  reduces  to 


u"kb||es“k  Vv(t^) 


+  •  •  •  l  +  •  • 


/Vkb\ 
;l  b-°\  Tkb  I 


+  . . .  I  +  ... 


(39b) 


and,  hence,  the  temperature  matches  if 


a  -  r1/(1+w) 

k  Tb 


/  i  \  -l/(l+w) 

l08(t)J  ^  1 


(40a) 


V* 


T—  .  , 


•  •  ^ 
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A,  0,  _  ■  0.  m  -  ^ 


k  k,®  b,o 


4’ 


(l  +  2c> 

Vn 

\  2  > 

s.  (l  +  1  *  *  ®b  ) 

b  \  a  b,w/ 

(A  Ob  ) 


The  PVL-BVL  longitudinal  velocity  matching  requires  that 


Vk  (skb-  -  ■••)  +  -}  '  k(v-  +  -) 


limkb,,<1  + 


.  1  +  .  •  •  j 

(41a) 


-  0 


i.e.,  that 


Umkb  II1  +  0,k 


Uk  (^) 

‘  Wb' 


l/(l+ui) 


+  .  .  . 


-  »„  — 

b'°  \  Tkb  / 


0  (41b) 


Thus,  the  longitudinal  velocity  matching  reconfirms  the  requirement 

*  -T1' 

k  b 


that  a,  m  Tl  (^+UJ),  and  yields  the  additional  requirement  that 


W.co  "  "Ub,o 
3  08 


m '  x  •  f  “»  <-2tk>  dtk  -  •  [(H^X1  +  v)"  ‘  • }] 

(4  2) 


The  PVL-BVL  radial  velocity  matching  reconfirms  the  requirement  that 
m  ■  3/4. 


An  important  fact,  implicit  in  the  results  of  Eqs.  (4Ga-b)  and  (42), 
is  that,  although,  to  the  order  of  approximation  considered  for  the  flow 
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quantities,  the  PVL-BVL  matching  requires  that  m  ■  3/4,  and,  hence,  the 
"effective  body"  shape  is  given  by 


This  matching  does  not  require  the  specification  of  the  (actual)  body 
shape,  and  so  the  equation  for  this  body  remains 


r  *  6^F^(x),  F^(x)  “  arbitrary  fnc(x) 


Therefore,  in  the  present  first-order  analysis,  the  flow  in  the  PVL  is 
independent  of  the  exact  body  shape.  This  point  is  not  made  explicitly 
in  the  first-order  analyses  of  Stewartson,  Solomon,  and  Ellinwood  and 
Mirels . 

Further,  since,  from  Section  IV, 


where 


1+U) 

Qf. 


a 


k 


1  /  (  1"Hd) 
Tb 


it  follows  that  the  parameters  and  are  given  (in  terms  of  the 
measure  of  the  actual  body  thickness  ratio,  6^,  etc.)  by 


(43a) 


1 
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(43b) 

In  addition,  the  inequalities  required  in  the  formulation  of  Eqs .  (28a -d), 


make  it  necessary  that  satisfy 

l0g(^)|  *  Tb/2  ^  1  (W) 

From  Eqs.  (A3)  and  (44),  6^,  itself,  must  satisfy 


VL<<: 


6b<<: 


« 


(45) 


which  requires  that  the  actual  body  must  be  very  slender  for  the  theory 
presented  to  be  valid.  For  example,  if  ■  10  then  6^  «  10  \ 

It  is  found  that  the  pressure,  friction,  and  heat  transfer  coef¬ 
ficients  are,  respectively, 


C 

P 


s.u. 
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(A6b) 


(A6c) 


l 


The  friction  drag  coefficient  is 


rxi  /^Al 

Jo  "i  H 

_ _ : — ILLh - 


3»u« 


(47a) 


The  pressure  drag  coefficient  is 


1 


-29- 


Xp 


/?  WA  (^rK 


P»ul  Cr^ 


C  V 
p,o  x 


lo8  I 


■1/2 


l^x'3'2 

[Fb(3T)> 

X 

2 

dx 

*  J  0 

L 

+  ... 


(A  7b) 


Thus,  the  principal  contribution  to  the  drag  coefficient  is  from 
friction,  since 


,  /^K1  ,v 

CD,£  \ 2c£,„/v 


x  !(v^ 


II  - 


lo8  |  ~ 2 

sb/ 


1/2 


ix-3/2f  7l/2 

*  n 


2 

*** 

X 

dx 

+  . . 


where 


—  «  1 
V 


In  terms  of  the  parameters  introduced  in  this  Memorandum,  the 

interaction  parameter  of  Ellinwood  and  Mirels,  x,  is  given  by 

(EM)  7 


‘(EM) 


x 

,2 


fkOO 


log  A 


(EM) 


lo8  |  ~~Z  I  +  ... 


(48) 


From  Eq.  (48),  then,  it  follows  that  the  expressions  for  C  C  C 

P*  f*  q’ 

CD,f  given  in  this  Memorandum  are  equivalent  to  those  found  by 
Ellinwood  and  Mirels  for  these  coefficients.  (Further,  since,  in  turn. 
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Tb  "  2  ^log  A(EM)]  +  ■*’ 

it  also  follows  that  the  expansion  parameter  of  the  PVL,  a^,  is  line¬ 
arly  proportional  to  the  expansion  parameter  of  Ellinwood  and  Mirels, 

e(EM)° 

When  the  (actual)  body  is  a  cone,  the  body  studied  by  Stewartson, 


Fb(x) 


x 


cone 


1 


and  the  pressure,  friction,  heat  transfer,  and  drag  coefficients  are 


The  foregoing  expressions  for 
by  Stewartson. 


and  are  equivalent  to  those  found 
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Appendix 

THE  VISCOUS  TRANSITION  LAYER 


To  span  the  distance  between  the  ISL  and  PVL,  a  viscous  transition 
layer  (TL)  is  introduced,  whose  coordinates  are  taken  to  be  and  rt, 
defined  by 


where  ia  a  parameter  much  less  than  unity  to  be  determined.  The 
expansions  for  the  flow  quantities  in  this  layer  are  taken  to  be 


u  -  1  +  ut  +  ... 

v  *  6{|  Nh  \lrt  +  Vt  +  •••} 


T 


S  k  t 


Tt  + 


P 


2  2 

ITS  Pt  + 


(A-2a ) 


The  behavior  of  u  and  T  in  the  ISL  ?VL  requires 


‘2  «  %  -  ■„ 


Tl/(1+U))  <<  t  (A- 2b) 


Then,  the  equations  of  motion  for  this  TL  reduce  to 
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dpt 

aT  "  °*  i-e“  pt  "  pt(xt>  " 


\Tt/  3xt\Tt/  4xt  Tt 


-  0 


(A-3) 


n 

To  find  a  self-similar  solution  for  T^,  take 


Tt  -  vt  et<ct> 


(»-4a) 


where 


V 


^t 


Vt 


(A-4b) 


where  a  is  a  constant  to  be  determined,  and  b  B  b(a)  *  3/A  +  [(1  +  u>)/2]a. 
Then,  the  continuity  and  energy  equations  combine  to  yield 


(A-5) 


For  the  sake  of  brevity,  self -similar  solutions  for  u£  and  are 
not  considered  in  this  presentation.  It  can  be  seen,  however,  from 
Eq.  (A-3),  that,  once  the  self -similar  solution  for  Tt  has  been  de¬ 
termined,  such  solutions  follow  directly. 


t 


1 
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* 

(This  equation  is  invariant  under  the  transformation  £  -«  £t  +  £t , 

* 

where  £  **  const  . ) 

Considering  that  the  TL  temperature  solution  should  approach  the 
PVL  temperature  solution  for  £  -*-“,  it  is  found  that  one  asymptotic 

solution  for  the  TL  temperature,  based  upon  a  balance  of  the  inviscid 
and  viscous  terms  of  Eq .  (A-5),  is 
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Considering  that  the  TL  temperature  solution  should  approach  the 
ISL  temperature  solution  for  £t  -♦  +  it  is  found  that  one  asymptotic 
solution  for  the  TL  temperature,  based  upon  the  domination  of  the  in¬ 
viscid  terms  in  Eq .  (A-5),  is 
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The  temperature  solutions  of  the  PVL  and  the  TL  are  matched 


through  the  introduction  of  the  intermediate  limit,  lim^,  defined  by 
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In  this  limit, 
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The  PVL-TL  temperature  matching  requires  that 
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From  inspection  of  Eqs.  (19),  (23),  and  (A-6),  it  can  be  seen  that 
the  matching  follows  directly. 

The  temperature  solutions  of  the  ISL  and  the  TL  are  matched 
through  the  introduction  of  the  intermediate  limit,  lim^t,  defined  by 
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In  this  limit, 
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The  ISL-TL  temperature  matching  requires  that 
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From  Eqs .  (14)  and  (A-7),  it  can  be  shown  that  the  matching  requires 
that 
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